Heat fluctuations in an out of equilibrium bath 
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We measure the energy fluctuations of a Brownian particle confined by an optical trap in an aging 
gelatin after a very fast quench (less than 1 ms) . The strong nonequilibrium fiuctuations due to the 
assemblage of the gel, are interpreted, within the framework of fluctuation theorem, as a heat flux 
from the particle towards the bath. We derive, from a simple model, an analytical expression of the 
heat probability distribution, which flts the experimental data and satisfles a fluctuation relation 
similar to that of a system in contact with two baths at different temperatures. 
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The heat flux between two reservoirs at different tem- 
peratures is an important and useful example of an out 
of equilibrium process. In small systems this heat flux is 
a strongly fluctuating quantity and the probability dis- 
tribution of these fluctuations has been recently widely 
studied, within the context of fluctuation theorems [1|. 
These studies have been mainly devoted to the steady 
state, that is when the temperatures Ta and Tb of the 
two reservoirs, A and 5, are kept constant. In such a case 
the probability distribution P{Qr) of exchanging with 
the reservoir A the heat Qr in a time r, is related to the 
that of exchanging the quantity — Qr according: 



In 



P{Qr) 
P{-Qr) 



Af3 Qr 



(1) 



where A/3 = (l/Tg — 1/Ta)//c^, /cb is the Boltzmann 
constant and kB^P Qr can be easily identified as the 
entropy production during the time r [ll-Q. This equa- 
tion has been derived for several theoretical models [H, Q 
in the stationary regime. However the non stationary 
case, although very useful for applications, has been stud- 
ied only in some specific models [5;] of systems relaxing 
towards equilibrium. Thus one may wonder whether a 
relation like Eq. ([T]) may still hold, how it is eventually 
modified and what kind of information on the system can 
be obtained [sj. 

These important questions have never been analyzed 
in any experiment. Thus the purpose of this letter is 
to give new insight to this problem, by measuring the 
energy fluctuations of a Brownian particle used as a probe 
inside a gelatin relaxing towards its solid-like state (gel), 
after a very fast quench, from above to below the gelation 
temperature Tg^i . The main result of our investigation is 
that these fluctuations can be interpreted as a heat flux 
from the particle towards the bath. The measured P{Qr) 
satisfies an equation formally equivalent to Eq. ([T]), but 
in this case A/3 is a decreasing function of time. The 
P{Qt) can be fitted by an analytical expression that we 
derive from a Langevin equation for a Brownian particle 
coupled with an out of equilibrium bath. 

The experiment has been performed using gelatin, a 




FIG. 1. (a) Schematic representation of the experimental 
setup to perform a local quench in a sol droplet around a 
trapped particle in the gel bulk, (b) Time evolution of the 
viscous drag coefficient 70 of the particle and the correlation 
time To of the gelatin droplet measured after the quench at 
/ = 5 Hz . Inset: kro/'jo as a function of time. 



thermoreversible gel obtained from denatured collagen. 
Above Tgei an aqueous gelatin solution is in a liquid vis- 
cous phase (sol), whereas below Tg^i the formation of 
a network of cross-linked filaments leads to an elastic 
sofid-like phase (gel) Q. In this gel phase the gelatin 
viscoelastic properties slowly evolve toward equilibrium 
and share some common phenomenological features with 
glassy dynamics [^-[sl • We are interested in this transient 
out-of-equilibrium regime, that we use to study the fluc- 
tuations of the energy fluxes from and to the heat bath 
in the nonstationary case. A similar problem has been 
theoretically analyzed for the first time in Ref. [5] for a 
model of aging spin glasses. It has been found that a re- 
lation like Eq. ([T]) can be applied to a relaxing system to 
obtain quantitative informations on the heat exchanges 
with the bath. We show in this letter that this approach, 
exploited only once in real experiments [10], can be in- 
deed very useful for understanding the properties of a 
Brownian particle in an out-of-equilibrium bath. 

In the present experiment, an aqueous gelatin solution 
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FIG. 2. (a) Time evolution of x after a quench. (Jx{t)^ is 
computed over = 0.1 s and over 60 independent quenches, 
(b) Time evolution of ax{t)^ (normalized hy ksT /k) after the 
quenches performed in gelatin (o) and glycerol (dashed line). 
Inset: Probability density of x at t = 0.5 s and 50 s for the 
quench in gelatin. The solid lines are Gaussian fits. 



(type-B pig skin) at a concentration of 10%wt is pre- 
pared following the usual protocol P. For this sample 
Tgei = 29° C. This solution fills a transparent cell the 
temperature of which is controlled by a Peltier element 
[Fig-^a)], at To = 26±0.05°C < Tgei. Thus the solution 
inside the cell is in the solid-like phase. A silica bead, of 
radius r = 1 /im, is kept inside the gelatin in the focal po- 
sition of a tightly focused laser beam (A = 980 nm) at a 
power of 20 mW. At this power the laser produces on the 
particle an elastic force of stiffness /c = 2.9 pN//im. Be- 
cause of light absorption, the temperature of the trapped 
particle is T = 27° C, which is still smaller than Tg^i. 
Therefore the bead is inside the solid gel in the beam 
focus at a distance h = 2b /im from the cell wall, see 
Figdja) [ll^. Starting from this condition, the laser 
power is increased to 200 mW and the local tempera- 
ture around the focus rises to 38° C > Tgei. As a result 
the gel melts and a liquid droplet of radius a = 5 /im, 
is formed around the trapped bead inside the the solid 
gel bulk, as sketched in Fig. [IJa). After 180 s, the laser 
power is suddenly decreased again to 20 mW so that the 
temperature is homogenized by heat diffusion into the 
bulk in less than 1 ms resulting in a very efficient quench 
of the droplet to the final temperature T < Tgei. At T 
the liquid inside the droplet solidifies in about 1 hour 



and the particle, trapped in the center of the drop by 
the focused beam, is a probe of this relaxation dynamics. 
The quenching procedure is repeated 60 times in order 
to perform the proper ensemble averages. 

Immediately after the quench we record the time evo- 
lution of the X position [see Fig. [U^a)] of the trapped 
particle measured by a position sensitive detector whose 
output is sampled at 8 kHz and acquired by a computer. 
The resolution of the measurement of x is about 1 nm 
lol . lij . In order to characterize the particle dynamics 
we measure, using active microrheology [lo|, the time 
evolution of the viscous drag coefficient 70 of the particle 
and the largest correlation time tq of the fluid. This is 
done by measuring the response of the bead at a time- 
dependent sinusoidal force F of amplitude 87 fN and fre- 
quency / applied to the bead. The force F = kxo is 
obtained through the modulation of the beam focus po- 
sition xq. The results for 70 and /ctq, measured at / = 5 
Hz, are shown in Fig. [TJb). First, for t < 200 s after the 
quench there is a transient regime where the droplet is 
purely viscous, tq — 0, whereas 70 increases in time. In 
this regime 70 and tq do not depend on /. For t > 200 s 
the liquid gelatine inside the drop has a behavior similar 
to that observed in macroscopic samples 0, , i. e. the 
liquid drop is actually undergoing gelation. We will study 
the nonequilibrium statistical properties of the bead dy- 
namics in the very first 200 s after the quench where the 
liquid gelatin inside the drop is mainly viscous and the 
elasticity is negligible with respect to /c, as shown in the 
inset of[T](b), where we plot /cro/70 as a function of time. 

We begin by analyzing the variance (Jx{t)'^ of x at time 
t after the quench. (Jx{t)'^ is computed over 60 indepen- 
dent quenches and over a short time window 5t = 0.1 s 
around each value of t in order to improve the statistics, 
as depicted in Fig. [2ja). The time evolution of (Jx{t)'^ is 
plotted in Fig. [Sfb). At the beginning, (Jx{t)'^ is almost 
three times the equipartition value ksT/k that would be 
obtained at equilibrium. This shows the presence of a 
stochastic force on the particle due to the transient for- 
mation of the gel network. This force weakens compared 
to the thermal fluctuations becoming negligible at ~ 20 
s so that cFxitY slowly decreases in time, reaching the 
equilibrium value for t > 20 s. This relaxation timescale 
is two orders of magnitude larger than the initial vis- 
cous relaxation time of the particle: r/e = 70//^ = 65 ms. 
Finally for t > 200 s, cFx(tY starts again to decrease be- 
cause of the appearance of a strong elastic component of 
the gel confirming the direct measure of 70 , tq , shown in 
Fig. mb), and justifying that for t < 200 s the gelatin 
elasticity is negligible. During this relaxation process x 
remains Gaussian as shown in the inset of Fig. EJb). 

In Fig. 2(b) we also plot the time evolution of cFx{tY 
measured, after the same quenching procedure, in a New- 
tonian fluid (glycerol 60%wt in water) with the same vis- 
cosity of the initial sol phase of gelatin. In this case, the 
particle dynamics must settle into an equilibrium state in 
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FIG. 3. (a) Probability density Pt{qr) of the normalized heat 
qt,T for T = 30 s at different times t after the quench. The 
theoretical lines are computed using Eq. ([5]). Inset: Pt{qr) 
computed, at the same t and r during a quench performed in 
glycerol. The solid line corresponds to the theoretical equi- 
librium profile, (b) Asymmetry function of the Ptiqr) shown 
in (a). The straight lines are obtained using Eq. (|6|). 



a time ^ Tk after the quench [l3|. Indeed in Fig.[2fb) we 
see that, in glycerol, (Jx{t)'^ = ksT/k for all t within the 
experimental accuracy. This confirms that no experimen- 
tal artifact is present and that the observed dependence 
of cFxitY gelatin is a real non-equilibrium effect due to 
the sol-gel transition. 

We now focus on the energy fluctuations of the particle 
inside the droplet for t < 200 s, i.e. when kro/jo ^ 1- 
Thus the energy of the particle is simply U {t) = kx{t)'^ /2. 
As there is no external force acting to the particle, the 
energy fluctuation AUt^r during the time r is equal to the 
heat Qt^T exchanged between the particle and the bath 
[l^ . specifically 



(2) 



(Qt,r) = (^/2)[cra;(t + T)2 -cr^(t)2] < 0, which reveals the 
existence of a mean heat flux from the particle to the sur- 
roundings over the timescale r because of the relaxation 
of (Jx. The maximum value \{Q^^^)\ ^ ksT takes place at 
t = s and for 20 s < r < 200 s. Non- negligible values 
of the mean heat compared to ksT persist for several 
seconds after the quench. Nevertheless, as t increases, 
I (Qt,r) I decreases becoming negligible and experimentally 
undetectable for t > 20 s . 

The probability density function Pt{qr) of the nor- 
malized heat qt^T = Qt.r/ {kBT) is computed over the 
60 quenches and over a short time window 5t = 0.1 s 
around each t and t + r, as sketched in Fig [2fa). We 
focus on a large value of r in order to probe timescales 
comparable to the relaxation of the nonthermal fluctu- 
ations. Fig. [3fa) shows Pt(^r) at different times t after 
the quench for r = 30 s. Ptiqr) is highly non- Gaussian 
with a spike at = and slowly decaying tails for all 
the values of t. Immediately after the quench, Pt{qr) is 
strongly asymmetric with a long tail occurring at nega- 
tive fluctuations. As t increases this asymmetry decreases 
and Ptiqr) becomes symmetric at t > 20 s.Once again, 
we check that the long-lived asymmetry occurs because 
of the intricate nonequilibrium nature of the bath. In the 
inset of Fig.[3fa) we plot Pt{qr) with r = 30 s for the local 
quenches performed in glycerol. Pt{qr) quickly converges 
to the equilibrium profile and it is always symmetric with 
respect to g^- = 0. 

As in Eq.([T]), the asymmetry function pt{qr) = 
ln[Pt{qr)/Pt{—qr))] is commonly used to measure the 
asymmetry of P{qr) between the positive and the neg- 
ative values of the fluctuations. The function pt{qr)^ 
computed from the Pt{qr) shown in Fig. [31(a), is plot- 
ted in Fig[3l(b). It is a linear function of its argument 
q^: pt{q^) = —A/St^rqr- The slope A/3t,r decreases as t 
increases approaching the symmetric value APt,r = as 
\{Qt,r) \ ^ ksT. This linear relation, except for the time 
dependent A/3t,r^ is formally similar to Eq.([T]), and it is 
the first experimental evidence of the phenomenon theo- 
retically obtained for a relaxing spin glasses in Ref. (Hf. 
For comparison we also plot pt{qr) for the quench in glyc- 
erol at t = and r = 30 s. In this case the heat exchange 
process is always symmetric, stressing that this is due to 
the nonequilibrium nature of the bath. 

In absence of a theory for our experimental results we 
model the nonequilibrium dynamics of the particle by an 
over damped Langevin equation for x: 



70 xt 



-kxt + 



for t < 200 s 



(3) 



The mean heat transferred during [t, t + r] is 



where —kxt is the harmonic force exerted by the optical 
trap, and ("t is a random noise representing the inter- 
action of the particle with the out of equilibrium bath, 
i.e. the drop undergoing gelation. Because of the out- 
of-equilibrium state, the statistical properties of C are 
unknown, and it is in general a nonstationary and corre- 
lated process. Multiplying Eq. (|3j) by Xt and integrating 
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over the time interval [t, t + r] one obtains the energy 
balance Eq. (j2j) jl^, where Qt^r is 



Q 
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rt+T nt+T 

J 70 ±5 ds + y ( 



(4) 



which is the sum of the viscous dissipation plus the the 
heat injected by the bath. Notice that Qt^r cannot be es- 
timated directly as ( is unknown. The only experimental 
way to measure Qt^r is via Eq. ([2]). 

The asymmetry of PtiQr) can be directly linked to the 
nonstationarity of the aging bath through the quantity 
Using Eq. (|2]) and the experimental fact that x is 
Gaussian [Fig. [2fb)], the analytical expression of PtiQr) 
for large r can be computed 



PtiQr) = ^^i^o(^t,rkr|)exp ( -^^^^^^Qr } , 



(5) 



where Kq is the zeroth-order modified Bessel function of 
the second kind, 



^t,T ~ /j- I \ ^ ^^^^ 



k(Jx(t)(Jx(t^T) 



and Bt^r = At^r^l^^lr/^ 

In Eq. ([5j) the asymmetry of the density is completely de- 
termined by the parameter A^^r in the exponential. At 
equilibrium At,r = 0, At^r = Bt.r = 1 regardless of t and 
r, so that one recovers the symmetric equilibrium profile 
Pt{qr) = Ko{\qr\)/7r with (^t,r) = [13. In Fig.lS^a) for 
each experimental Pt{qr) we plot the theoretical predic- 
tion given by the analytical formula ^ using the respec- 
tive experimental values of cix shown in Fig. EJb). The 
excellent agreement confirms that the Langevin model 
([3]) is suitable to describe the particle dynamics and the 
heat exchange with the gelatin bath after the quench. 

From Eq. (|5]) one obtains the explicit expression for 
the asymmetry function ptiQr) = —^Pt^r^r and Al3t,T 



k 



(6) 



Hence, the linearity of ptiQr) is analytically satisfied 
for all the values of the heat fiuctuations and for all t 
even when PiQr) is strongly non-Gaussian. In Fig. [3jb) 
we plot the straight lines with the slope AfSt^r given 
by Eq. (|6j) and computed using the experimental val- 
ues of The good agreement with the experimental 
data shows that Eq. (|5j) verifies a fiuctuation relation, as 
Eq. (dl). 

Eq.(j6]) gains a very intuitive interpretation if one intro- 
duces an equipartition-like relation for the particle mo- 
tion for < t < 60 s: ksTeffit) = kcjxitf . Here Tg// 
is the effective temperature perceived by the particle due 



to its coupling with the nonequilibrium gelatin environ- 
ment. In this way the parameter APt,T can be written 
conveniently as Aft,^ = [l/Te//(t + r) - l/Te//(t)] T, 
which is formally equivalent to that of Eq. ([T]). Hence 
ASr = —kBAPt,TQt,T can be naturally identified as the 
entropy produced by the breakdown of the time-reversal 
symmetry due to the the effective temperature imbal- 
ance at two different times after the quench. As the 
gelatin droplet ages AS'^- slows down and the particle 
exhibits an equilibrium-like dynamics for the experimen- 
tal timescales. We point out that unlike Eq. ([1]) derived 
in Refs. [H, for nonequilibrium steady states, Eq. ([6]) 
holds for a nonstationary regime created by the nonequi- 
librium bath. 

In conclusion, we have experimentally studied the fiuc- 
tuations of the heat exchanged between a trapped Brow- 
nian particle and a non-stationary bath, i.e. an aging 
gelatin after a very fast quench. We have shown that 
the distribution of the heat satisfies a fiuctuation rela- 
tion even when the bath is in a non-stationary state. A 
Langevin model justifies the observation. The analogy of 
our results with those obtained for spin glasses suggests 
that this fiuctuation relation may appear as a very robust 
symmetry property of heat exchange processes in other 
kinds of relaxing systems. 
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